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We analyze pp total cross section data from accelerator to cosmic-ray energies and show that the 
experimental information presently available indicate two different scenarios at the highest energies. 
One of them accommodates the predictions of the majority of models and the other one shows 
agreement with a multiple diffraction model recently developed and improved. We argue that 
these results may bring a kind of warning against some possible precipitated attempts to constraint 

- - . cosmic-ray estimations from analyzes of accelerator data. 
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Some years ago, a successfull description of pp elastic scattering data in the interval 13.8 < y 7 ! < 62.5 GeV was 
£>^' obtained through a multiple diffraction model |Q, hereafter referred as MDM. Extrapolations of the predicted af£ t 
showed agreement with the highest quoted values of the cosmic-ray experiments. In particular, the model predicts 
a tot(\/~s — 16 TeV) = 147 mb (without error estimated) |Q. To our knowledge this is the only phenomenological 
approach that indicates that the cross section may rise faster with the energy than generally expected from "standard" 
i model predictions |^-||] and some fit results Q . 

Recently, Perez-Peraza et al. introduced some improvements in the model predictions by determining confident 
error bands, through a forecasting regression analysis [7j]. According to the authors, the extrapolations present so large 
error bands that the predictions show agreement with all the cosmic-ray results, including those that are incompatible 
, among themselves. Moreover, the authors show that taking account of the pp data in the interval 13.8 < y/s < 62.5 
GeV together with pp data at 546 GeV and 1.8 TeV, the predictions resulting from a new fit of the model free 
parameters present a very narrow error band and are compatible with the smallest quoted values of the cosmic-ray 
experiments. Since this seems to be in agreement with "standard" models and fit results the authors conclude 

that "extrapolations from accelerator data should be used to constraint cosmic-ray estimations" 0, which was also 
suggested in Q and done in practice in ||] . 

However some aspects of the approach must be commented. Firstly, reading from the plot in Fig. 2 of this 
reference we can infer af^ t (^/s = 16 TeV) ~ 147 ± 37 mb. Although this is a large error band (~ 25%) the results still 
favor (at least qualitatively) the highest quoted values of the cosmic-ray experiments, as shown in the same Figure 
(this is not what the authors claimed). Secondly, in the context of a model approach, simultaneous analysis of pp and 
i-C • pp scattering should take analyticity and crossing into account, what was not done in and also in 0. At last, the 
assumption that pp and pp interactions are the same even at 546 GeV is a strong ad hoc hypothesis and since pp and 
pp interactions are different in the ISR energy region, one should expect the use of af^ t and not of£ t in this interval. 

In order to attempt to bring some insight on this issue, we performed a simple model-independent analysis of the 
same pp accelerator data used in [Q and but including cosmic-ray information presently available. As reviewed 
in some detail in jjj, at cosmic-ray energies (-y/s : 6 — 40 TeV) we can identify two distinct set of estimations, one 
represented by the results of the Fly's Eye Collaboration |jj together with those by the Akeno Collaboration flio| ], 
and other by the results of Gaisser, Sukhatme, Yodh (GSY) |ll|] with those by Nikolaev [Jil l . All these data and 
estimations are displayed in Fig. 1 (Akeno results were read from the graph of Fig. 4 in Ref. JlOfl). 

The discrepancies between these two sets of cosmic-ray information comes from both, the uncertainties in the 
determination of the production cross section <7p™ a d ir , and the uncertainties in the connection between o^° air and the 
pp cross section, as discussed in detail by Engel, Gaisser, Lipari and Stanev all the cosmic-ray estimations are 
strongly model dependent. One important point is the fact that the results by Nikolaev and GSY have never been 
directly critized in the literature. 

For these reasons, we shall investigate the behavior of the total cross section by taking into account the discrepancies 
that characterize the cosmic ray information. Also, since the base of our discussion is the predictions of the MDM, 
we shall consider the same set of pp accelerator data used in references ]l| and {^]. In order to accomplish that, we 
select three ensembles of data and experimental information with the following notation: 

• Ensemble I: pp accelerator data [[lJJ , 

• Ensemble II: pp accelerator data + Akeno + Fly's Eye (l4||^,[I(| , 

• Ensemble III: pp accelerator data + Nikolaev + GSY [pUlJO]. 
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FIG. 1. Proton-proton total cross sections: accelerator data (l4| and cosmic-ray estimations |] |l^| 

Different analytical parametrizations may fit all these ensembles with quite good statistical results, as shown, for 
example, in Jlql . In this communication we shall choose the simplest form, which corresponds to the same analytical 
function used to fit the free parameters, depending on the energy, in the MDM. As explained in JlJ this concerns a 
polynomial of second degree in Ins. In terms of of^ this choice satisfies the Froissart-Martin bound, 



<rZ = A + B\n- + C 



s 

•so 



ln± 

so 



(1) 



where A, B, C are free parameters and so = 1 GeV 2 . 

It should be stressed that in JlJ the fits involved only the differential cross section and the p parameter data (input 
information). The total cross section was predicted from parametrizations of essentially two free parameters depending 
on the energy, one associated with the hadronic form factor and the other with the average elementary scattering 
amplitude. Since these parameters enter in the eikonal (in the momentum transfer space), the final result for af^ t 
should not be expected to be the same as that obtained with Eq. (1). In addition, it should be noted that other 
parametrizations including terms like (Ins) 13 and Rs^ 1 / 2 (D and R free parameters) lead to a faster rise of a^ t (s) 
than that obtained through Eq. (1), as shown in 15|. In this sense the above choice clearly constraint the rise of crf^ 



and shall be considered as a lower limit (allowing some connection with the parametrizations of the MDM) . 

The fits have been performed using the CERN-MINUIT routine The corresponding values of the free pa- 

rameters and statistical information about each fit are displayed in Table 1. Errors from the free parameters were 
propagated to ofo t (s) through the usual formulas, taking account of both variances and covariances |P7j ]. The results 
with the corresponding error bands (± one standard deviation) for ensembles I, II and III are shown in Figs. 2, 3, 
and 4, respectively. 



Table 1. Values of the parameters A, B and C in Equation (1) (all in mb) and the x 2 per degree of 

freedom in each fit to Ensembles I, II and III. 



Ensemble: 



I 



II 



III 



A 
B 

C 

x 2 

d.o.f. 
X 2 /d.o.f. 



47.32 
-3.809 
0.4042 
8.654 

4 
2.16 



45.78 
-3.315 
0.3654 
11.75 
11 
1.07 



49.27 
-4.435 
0.4534 
15.04 
11 
1.37 
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From Fig. 2, the results with the ensemble I lie between the two sets of discrepant cosmic-ray estimations. It 
corresponds to extrapolations from a "pure" accelerator information in the interval investigated. 

On the other hand, if one or the other cosmic-ray sets are included the situation changes, as shown in Figs. 3 and 
4. Taking into account the error bands (one standard deviation), the two fits become discrepant above ~ 1 TeV and 
they are certainly incompatible above 5 TeV. We understand this result as a statistical evidence for the possibility of 
two distinct scenarios above yfs ~ 5 TeV. 

In particular, from each ensemble, the interpolated results at y/s = 14 TeV are af% t = 125 ±7 mb (I), afj t = 119 ±5 
mb (II), and of^t — 133 ± 6 mb (III). In Ref. jlj| four different parametrizations were used to fit Ensembles II 
and III and the two scenarios are also inferred. For example, taking into account the 4 parametrizations (without 
error estimations) we can infer from Ensemble II the value a^ t (14 TeV) = 113 ± 5 mb and from Ensemble III, 
ofo t (14 TeV) = 140 ± 7 mb. At lower values of energy, namely \fs in the region 10 — 100 GeV, there is no significant 
distinction between the four parametrizations. 
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FIG. 2. Parametrization to Ensemble I: pp accelerator data 
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FIG. 3. Parametrization to Ensemble II: pp accelerator data + Akeno + Fly's Eye |l4|,P,[lo[ . 
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Ensemble II clearly accommodates the predictions from the majority of models and in general is considered as "the 
cosmic-ray results" . On the other hand, if the estimations in Ensemble III are the correct ones, the experimental 
information on a^ t presently available indicate that the pp cross section rises faster with the energy than generally 
expected or assumed and, in this case, the rise is well described by the MDM. At this point it should be stressed that, 
due to its phenomcnological basis, this model can not provide a "microscopic" explanation for this phenomena (only, 
perhaps, some naive geometrical arguments) and this puts limitations in the understanding of a "fast rising" . 

However, it should be remembered that several cosmic-ray experiments point out to the possibility of new phenomena 
in pp collisions at yfs ~ 5 TeV and this has never been investigated by accelerators. Moreover, recently, more 
confident analyzes stress the particular features of Centauro events, when compared with normal events observed in 
emulsion chamber experiments | fl9|| . Since the bulk of the cosmic-ray showers concerns pp and not pp interactions, 
these "exotic" events may represent an additional contribution to the pp channel, resulting in a faster increase of the 
ofot i as claimed in po[ . 

Certainly, these puzzles shall be resolved by the CERN Large Hadron Collider (LHC) and even the BNL Relativistic 
Heavy Ion Collider (RHIC). But before this, we understand that the results presented here may serve as a kind of 
warning against possible precipitated attempts to constraint cosmic-ray estimations from analyzes of accelerator data. 

Since in this communication our central goal was to discuss the predictions of the MDM, the same set of pp 
accelerator data was used. An extend analysis with all experimental data and estimations available, including pp 
scattering, is in course. 
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